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SUMMARY 


The  compounding  technique  of  obtaining  stress  intensity  factors  is  applied 
to  cracks  at  the  edge  of  one  hole  in  a row  of  periodically  spaced  holes.  The 
cracks  (one  or  two)  lie  along  the  line  of  holes  and  perpendicular  to  the  uniform 
stress  which  is  applied  remote  from  the  holes.  This  configuration  models  cracks 
at  rivet  holes  in  airframe  structures  reinforced  with  riveted  stiffeners,  if 
there  is  little  in-plane  load-transfer  through  the  rivets.  A particular  feature 
of  the  solution  is  the  high  accuracy  of  the  results  at  short  crack  lengths  - a 
necessary  requirement  for  the  reliable  estimation  of  the  lifetimes  of  growing 
cracks . 
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1 INTRODUCTION 

In  engineering  materials  used  in  aircraft  structures  cracks  may  be  present 
either  at  manufacture  or  they  may  appear  during  the  service  life.  To  ensure 
safety  and  to  optimise  inspection  schedules,  it  is  necessary  to  know  both  the 
residual  strength  of  the  cracked  structure  and  the  rate  at  which  the  crack  grows 
when  subjected  to  service  loads.  Both  are  dependent  on  the  stress  intensity 
factor  which  governs  the  stress  field  at  the  tip  of  the  crack. 

1-3 

The  stress  intensity  factor  is  known  for  many  simple  configurations 
However,  in  practice,  cracks  occur  in  more  complex  configurations.  Most  methods 

4 

for  determining  stress  intensity  factors  , either  theoretically  or  experimentally, 
are  both  costly  and  time-consuming.  Hence  there  is  a need  for  a simple,  quick, 
and  easy-to-use  method  for  obtaining  stress  intensity  factors  for  practical 
configurations  from  the  known  solutions  for  simple  ones.  A compounding  method 
was  developed  by  Cartwright  and  Rooke^  in  which  boundary-boundary  interactions 
are  not  significant.  This  method  has  been  modified  to  include  cracks  in  stiffened 
structures^  and  further  modified  to  include  the  case  of  a crack  at  the  edge  of  a 
hole^.  In  this  latter  modification,  it  was  necessary  to  include  a contribution 
for  boundary-boundary  interaction. 

In  this  Report  the  compounding  method  is  applied  to  the  problem  of  cracks 
at  the  edge  of  a hole  which  is  in  a row  of  holes  in  a uniformly  stressed  sheet; 
the  stress  which  acts  in  a direction  perpendicular  to  the  row  of  holes  causes 
the  cracks  to  grow  along  the  line  of  holes.  This  configuration  is  representative 
of  some  airframe  components.  For  example,  cracks  may  develop  at  a rivet  hole  in 
a longitudinal  stiffener  in  a pressure  cabin,  in  which  the  major  in-plane  loading 
arises  from  the  hoop  stress,  with  little  in-plane  load  transfer  through  the 
rivet.  An  analagous  situation  may  arise  in  a wing  skin  in  which  a chordwise 
crack  may  develop  at  a rib/skin  fastener  under  wing  bending  fatigue  loads. 

In  order  to  apply  the  compounding  method  to  such  configurations  the  stress 

intensity  factors  are  required  for  the  following  ancillary  configurations;  an 

isolated  hole  with  one  or  two  cracks  (Case  1.3.3,  Ref  I)  and  a crack  near  an 

isolp^'ed  hole  (Case  1.3.5,  Ref  1).  The  boundary  interaction  term  is  estimated 

fro  I’esults  for  a cracked  hole  subjected  to  localized  forces  on  its 

pe:  ’ . The  compounding  method  gives  accurate  solutions  for  short  crack- 

lengtns  * the  stress  intensity  factor  tends  to  that  of  a crack  at  the  edge  of  a 

hole  with  a stress  oonaentvation  factor  K , which  is  determined  by  the  geometry 

9 10  ^ 

of  the  uncracked  structure  ' . The  importance  of  accuracy  for  short  cracks  lies 
in  the  fact  that  a major  portion  of  the  structures  useful  life  is  spent  when 
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the  crack  is  short.  Unless  reliable  predictions  of  strength  and  crack-growth 
rate  of  short  cracks  can  be  made,  safety  considerations  may  result  in  the  use  of 
over-conservative  assumptions  in  design  and  inspection  frequency,  with  consequent 
economic  penalties. 

2 THE  COMPOUNDING  METHOD 

The  method  of  compounding  is  described  in  Ref  5 where  the  resultant 
normalized  stress  intensity  factor  is  found  by  using  ancillary  configurations 

that  in  combination  approximate  to  the  one  under  consideration.  The  structural 
configurations  studied  in  this  Report  are  those  of  infinite  sheets  with  rows  of 
holes  (radius  R)  in  a line,  uniformly  stressed  remote  from  and  perpendicular 
to  the  line  of  holes  (see  Fig  1).  Two  configurations  are  considered;  the  first 
where  one  of  the  holes  is  symmetrically  cracked  and  the  second  where  the  hole 
has  only  one  crack.  In  both  cases  a is  the  distance  from  the  centre  of  the 
hole  to  the  tip  of  the  crack.  In  order  to  calculate  the  effect  of  other 
boundaries  on  a cracked  hole  we  have  to  substitute  for  the  cracked  hole  an 
equivalent  simple  crack,  which  has  the  same  normalized  stress  intensity  factor 
as  the  cracked  hole;  the  procedure  used  is  the  same  as  that  in  Ref  7.  In  order 
to  obtain  an  expression  for  the  equivalent  crack  length,  we  consider  the  stress 
intensity  factor  of  a simple  crack  of  length  2a'  in  an  infinite  sheet  under 
uniform  uniaxial  tensile  stress  a remote  from  the  crack, 

K ■ o/naT  (1) 

as  a definition  of  the  equivalent  crack.  This  must  be  the  same  as  that  for  an 
isolated  cracked  hole,  which  is  and  is  given  by, 

Kq  - QQo/ira  , (2) 

where  Qq  is  the  normalised  stress  intensity  factor  of  a cracked  hole  remote 
from  all  boundaries.  Comparison  of  equations  (I)  and  (2)  gives 

a’  - Q^a  . (3) 

This  means  that  for  the  ancillary  configurations  the  cracked  hole  is  replaced  by 

2 

a crack  of  length  2QQa  (see  Fig  2) . 
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vnien  the  cracked  hole  is  replaced  by  an  equivalent  crack  of  length  2a' 
the  distance  from  the  crack-tip  to  the  boundary  is  kept  constant,  ie  for  the 
right-hand  tip  (tip  A in  Fig  2) 


b'  - a' 
n 


b - a 
n 


(4) 


where  b^  is  the  distance  from  the  centre  of  the  equivalent  crack  to  the  nth 


(n  > 0)  boundary  on  the  right  and  b^  is  the  distance  from  the  centre  of  the 


cracked  hole  to  the  nth  boundary;  similarly  for  the  boundaries  on  the  left 
(n  < 0) 


b’  + a'  - 
n 


b + a 
n 


(5) 


The  basic  compounding  formula  for  a crack  in  the  presence  of  an  infinite 


array  of  boundaries,  is  given  by  , 


n>co 


K - K + 
r 


I 


+ K 


(6) 


n*-* 


where  R is  the  stress  intensity  factor  of  a crack  of  length  2a  in  an  infinite 
sheet  subjected,  remote  from  the  crack,  to  a uniform  uniaxial  tensile  stress  a 
and  R is  given  by 

a/it&  . (7) 


R - 


R is  the  stress  intensity  factor  of  the  crack  in  the  presence  of  the  nth 
n 


boundary  only  and  R is  the  contribution  due  to  boundary-boundary  interactions. 


Equation  (6)  can  be  written  in  terms  of  the  normalized  stress  intensity 
factor  Q(«  R/R)  as 

n»« 


1 + 


1)  ♦ Q, 


(8) 


n»  -» 


The  introduction  of  the  equivalent  crack-length  a'  into  the  ancillary 


configurations  modifies  equation  (8)  to 


y 


1 + 


i 


1) 


n*-» 


♦ Q. 


* 0 


(9) 


where  the  normalized  stress  intensity  factors  are  for  an  isolated  cracked 

hole  and  for  the  equivalent  crack  in  the  presence  of  the  nth  boundary  only. 

A general  expression  for  the  correction  has  been  derived^  using  an 

alternating  technique;  it  is  not,  however,  easy  to  calculate.  It  has  been 

suggested^  that  an  approximation  for  may  be  obtained  by  considering  the 

total  effect  to  be  that  of  two  point  forces  P applied  to  the  top  and  bottom  of 

the  cracked  hole.  The  magnitude  of  the  forces  P is  determined  so  as  to  give 

the  same  stress  oonoentvation  factor  at  the  edge  of  the  uncracked  hole  as 

that  which  occurs  due  to  all  the  other  boundaries  in  the  configuration.  These 

forces  P are  a function  of  the  individual  configuration  through  its  K which 

9 10  ^ 

is  available  for  many  common  configurations  * . Values  of  , for  a hole 

with  one  or  two  cracks  as  a function  of  the  crack  length  divided  by  the  hole 
radius,  are  given  in  Ref  7,  where 


2P 

irRo 


where  is  the  normalized  stress  intensity  factor  as  a* 

of  the  nth  boundary  only. 

The  complete  coiiq>ounding  formula  is  therefore; 


(10) 


0 in  the  presence 


^r  - ‘^O 


r » 1 

/ “ \ 

+ 

1 

+ 

*t-  (m;-  1)1 

L n-l 

\ / 

(11) 


2. 1 Application  to  a symmetrically  cracked  hole 

The  geometry  of  the  configuration  is  defined  by  the  hole  radius  R , the 
distance  a between  the  hole-centre  and  the  crack-tip,  and  the  spacing  b 
between  hole-centres. 

Graphs  of  and  for  a single  hole  with  two  equal-length  cracks 

are  given  as  a function  of  a/R  in  Ref  7.  Values  of  for  the  ancillary 

configuration  of  a crack  located  near  a hole  can  be  obtained  from 
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Case  1.3.5,  Ref  I.  is  given  as  a function  of  ^n^^  where 

is  the  distance  from  the  centre  of  the  equivalent  crack  to  the  edge  of  the  nth 
hole  in  the  ancillary  configuration  (d^  > b^'  > R).  The  dependence  of  ^'/d^ 
and  dg/R  on  n/R  nnd  b/R  are  given  in  the  Appendix.  The  first  term  in 
equation  (II)  can  be  evaluated  from  known  values  of  Qq  and  ; the  value  of 
(0^  - I)  was  negligible  for  |n|  > 2 for  all  values  of  a/R  and  b/R  . 

Values  of  and  are  required  in  order  to  calculate  the  boundary 

interaction  term  Q . The  values  of  K.  for  various  b/R  ratios  are  given^’*^ 
e v 

in  Table  I . 


Table  I 

Stress  oonoentration  factors  for  a hole 
in  a periodic  row  of  holes 


b/R 

3.0 

3.5 

4.0 

5.0 

10 

3.92 

3.44 

3.24 

3.10 

3.01 

The  values  of  M*  are  obtained  from 
n 


M'  ■ Liffl  {QJ*}  ; 
“ a*-K)  ^ 


(12) 


the  appropriate  values  of  n'/d^  and  <1q/R  nre  given  in  the  Appendix.  For  the 
configurations  considered  here  the  boundary  interaction  term  was  small 
(Qe  < 5X  Qj). 

Curves  of  K^/(ai^rR)  , ie  /a/R  , are  plotted  as  a function  of  ^ a/R  , 
in  Figs  3 and  4,  for  various  values  of  b/R.  Fig  3 contains  the  results  for 
short  cracks  I.O  4 a/R  4 I. 25  and  Fig  5 for  longer  cracks.  In  Fig  4 the 
stress  intensity  factor  tends  to  that  for  a line  crack  of  length  2a  located 
symsetrically  in  an  array  of  holes;  the  solution  for  this  configuration  may  be 
obtained  by  compounding  from  Case  1.3.7,  Ref  I.  Two  curves  are  shown  for  this 
approximation  (b/R  ■ 3 and  4)  in  Fig  4. 

2.2  Application  to  a hole  with  one  crack 

The  geometry  of  the  configuration  is  defined,  as  before,  by  the  hole 
radius  R , the  distance  a from  the  hole-centre  to  the  crack-tip  and  the 
spacing  b of  the  holes. 
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The  functions  Q.  and  Q /&  which  are  needed  to  calculate  Q from 

° 1 ® ^ . 
equation  (II)  are  available  , as  a function  of  a/R  . The  ancillary  configuration 

is  again  the  equivalent  crack  in  a row  of  holes  (Fig  1);  the  parameters 

and  which  are  required  to  obtain  (Case  1.3.5,  Ref  I)  are  given  in 

the  Appendix.  Again  (Q^  - 1)  was  negligible  for  |n|  > 2 . Since  and 

are  evaluated  for  zero  crack-length,  they  are  the  same  as  in  the  two-crack 

proolem. 

The  resultant  normalized  stress  intensity  factor  can  now  be  compounded 

from  equation  (M);  the  boundary  interaction  term  was  again  small  (Q^  < 5%  Q^.)* 
Curves  of  K^/(ai^i^)  , ie  »^a/R  , are  plotted  in  Figs  5 and  6 as  a function  of 

a/R  for  various  values  of  b/R  . Fig  5 contains  the  results  for  short  cracks 
1.0  4 a/R  j 1.25  and  Fig  6 for  longer  cracks.  For  the  longer  cracks  the  stress 
intensity  factor  will  tend  to  that  for  a line  crack  of  length  a R located 
eccentrically  in  an  array  of  holes;  the  solution  for  this  configuration  may  be 
obtained  by  compound'ng  from  Case  1.3.5,  Ref  1. 

3 DISCUSSION 

In  this  Report  the  compounding  technique  has  been  used  to  calculate  the 
stress  intensity  factor  for  a cracked  hole  in  a row  of  equally  spaced  holes.  The 
technique  is  not  limited  to  the  configurations  considered,  ie  one  crack  or  two 
equal-length  cracks  at  the  edge  of  one  of  the  holes.  The  more  general  problem 
of  a hole  with  two  unequal-length  cracks  in  a non-periodic  row  of  dissimilar  holes 
can  be  solved  using  the  same  ancillary  configurations^’^.  The  presence  of  the 
holes  near  the  cracked  hole  increases  the  stress  intensity  factor  over  that  of  an 
isolated  cracked  hole;  the  effect  is  negligible  if  the  holes  are  more  than  lOR 
apart,  and  increases  as  the  hole-spacing  decreases  (see  Figs  3 to  6).  In  all 
cases  the  effect  of  the  hole  nearest  to  the  tip  being  considered  is  about  ten 
times  greater  than  the  effect  of  the  next  nearest  hole  on  the  same  side  of  the 
tip;  the  effect  of  all  the  other  holes  on  the  same  side  of  the  tip  is  negligible 
(4  IZ).  The  holes  on  the  far  side  of  the  crack  from  the  tip  being  considered 
have  a smaller  effect  (20  to  25Z)  than  those  on  the  near  side. 

A particular  advantage  of  the  compounding  technique  for  cracked-hole 
problems  is  that  the  stress  intensity  factor  is  determined  accurately  for  short 
cracks.  This  accuracy  is  ensured  by  relating  the  boundary-boundary  interactions 
to  the  stress  oonoentration  factor  K in  the  uncracked  configuration.  For 
long  cracks  the  solutions  obtained  approach  those  obtained  by  replacing  the 
hole  plus  cr8ck(8)  by  a simple  crack  of  the  same  overall  length  as  the  hole  plus 
crack(s).  It  is  estimated  that  the  maximum  errors  in  the  stress  intensity 
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factors  are  4102  overall,  and  much  less  for  the  important  region  of  short  cracks; 
it  is  the  behaviour  of  short  cracks  which  dominates  some  damage^tolerant  design 
calculations  since  a large  proportion  of  the  useful  life  of  a cracked  component 
is  spent  while  the  crack  is  short. 

Although  the  boundary-boundary  interaction  term  was  included  in  order 

to  ensure  accuracy  for  short  cracks  the  contribution  to  the  stress  intensity 
factor  was  small  (^5Z)  for  the  configurations  considered.  A more  closely  spaced 
array  of  holes  would  result  in  larger  values  of  . 

Comparison  of  Figs  3 and  A with  Figs  5 and  6 show  that  for  the  same  values 
of  b/R  and  a/R  the  stress  intensity  factors  for  one  crack  are  always  less 
than  for  two  cracks;  for  short  cracks  (a/R  < 1.25)  the  difference  is  less  than  SZ. 

A CONCLUSIONS 

(1)  The  compounding  technique  can  be  applied  to  a cracked  hole  in  a row  of 
holes  to  give  stress  intensity  factors  for  a conmon  airframe  configuration. 

(2)  The  solutions  are  acceptable  for  most  engineering  applications;  they  are 
particularly  accurate  for  small  cracks,  a necessary  condition  for  the  use  of  the 
method  in  damage  tolerant  designs. 
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Appendix 


DETERMINATION  OF  PARAMETERS  FOR  ANCILLARY  CONFIGURATIONS 

(See  section  2) 


The  parameters  required  in  the  ancillary  configuration 


must  be  expressed  in  terms  of  a/R  and  b/R,  The  condition  that  the  distance 
from  the  crack-tip  to  the  boundary  should  be  the  same  in  the  ancillary  configura- 
tion as  in  the  original  configuration  led  to  equations  (4)  and  (5)  which  can  be 
written  in  terms  of  d the  distance  from  the  crack  centre  to  the  boundary  edge: 


and 


d' 

n 


d - a + a',  n>0 

n 


d’ 

n 


d +a-a',  n<0 

n 


(A-1) 


Since 


n 


|n|b  - R , 


(A-2) 


equation  (A-1)  becomes 


d' 

n 


~ - 1 + — (0^  - 1) 


n > 0 


and 


> 


d' 

n 


Mk  - I _ £ - n 

R^  ' R % * 


n < 0 . 


(A-3) 


Therefore,  since  a'  - Q^a  , 


a 

^0  R 


TTr 

n 


(A-4) 


For  the  special  case  a'  ■ 0 , required  in  the  calculation  of  M’  , the  parameter 


d^/R  have  the  following  values: 


nb 


~ - 2 , for  n > 0 , 
n K 


and  d'/R 
n 


nb 


, for  n < 0 . 


....  (A-5) 
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